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Abstract 

We analyze lepton asymmetry induced in the decay of right-handed neutrinos in 
a class of minimal left-right symmetric models. In these models, which assume low 
energy supersymmetry, the Dirac neutrino mass matrix is proportional to the charged 
lepton mass matrix. As a result, lepton asymmetry is calculable in terms of 9 param- 
eters, all measurable in low energy neutrino experiments. By solving the Boltzmann 
equations numerically we show that adequate baryon asymmetry is generated in these 
models. This however places significant constraints on the light neutrino parameters. 
We find tan 2 — m\/rri2 and #13 = (0.01 — 0.07) for the neutrino oscillation angles, 
and (3 ~ a + ir/2 for the Majorana phases. 
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1 Introduction 



The discovery of neutrino flavor oscillations in solar, atmospheric, and reactor neutrino ex- 
periments pQ may have a profound impact on our understanding of the dynamics of the 
early universe. This is because such oscillations are feasible only if the neutrinos have 
small (sub-eV) masses, most naturally explained by the seesaw mechanism [2]. This as- 
sumes the existence of super-heavy right-handed neutrinos AT f (one per lepton family) with 
masses of order (10 8 — 10 14 ) GeV. The light neutrino masses are obtained from the matrix 
M v ~ MpM^ 1 Md T where Md and Mr are respectively the Dirac and the heavy Majo- 
rana right-handed neutrino (r.h.n) mass matrices. The decay of the lightest right-handed 
neutrino Ni can generate naturally an excess of baryons over anti-baryons in the universe 
[3| consistent with cosmological observations. The baryon asymmetry parameter is an im- 
portant cosmological observable constrained by Big Bang Nucleosynthesis and determined 
recently with high precision by the WMAP experiment [3]: 



The decay of N\ can satisfy all three of the Sakharov conditions [5] needed for successful 
generation of T]b - it can occur out of thermal equilibrium, there is sufficient C and CP 
violation, and there is also baryon number violation. The last condition is met by combining 
lepton number violation in the Majorana masses of the right-handed neutrinos with B + 
L violating interactions of the Standard Model arising through the electroweak sphaleron 
processes 0. A compelling picture emerges, with the same mechanism explaining the small 
neutrino masses and the observed baryon asymmetry of the universe. t)b appears to be 
intimately connected to the observed neutrino masses and mixings. 

A more careful examination of the seesaw structure would reveal that, although there is an 
underlying connection, the light neutrino mass and mixing parameters cannot determine the 
cosmological baryon asymmetry, when the seesaw mechanism is implemented in the context 
of the Standard Model (SM) gauge symmetry. It is easy to see this as follows. Without 
loss of generality one can work in a basis where the charged lepton mass matrix and the 
heavy right-handed neutrino Majorana mass matrix Mr are diagonal with real eigenvalues. 
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The Dirac neutrino mass matrix would then be an arbitrary complex 3x3 matrix with 18 
parameters (9 magnitudes and 9 phases). Three of the phase parameters can be removed by 
field redefinitions of the left-handed lepton doublets and the right-handed charged lepton 
singlets. The neutrino sector will then have 18 (= 15 + 3) parameters. 9 combinations of 
these will determine the low energy observables (3 masses, 3 mixing angles and 3 phases), 
while the lepton asymmetry (and thus r] B ) would depend on all 18 parameters, leaving it 
arbitrary. 

In this paper we show that it is possible to quantitatively relate t]b to light neutrino mass 
and mixing parameters by implementing the seesaw mechanism in the context of a class of 
super symmetric left-right models 7j. We note that unlike in the SM where the right- 
handed neutrinos appear as rather ad hoc additions, in the left-right symmetric models they 
are more natural as gauge invariance requires their existence. Supersymmetry has the well- 
known merit of solving the gauge hierarchy problem. With the assumption of a minimal 
Higgs sector, it turns out that these models predict the relation for the Dirac neutrino mass 
matrix 

fm e \ 



M, 



D 



m M 



(2) 



\ m T J 

where c ~ m t /mb is determined from the quark sector, leaving only the Majorana mass 
matrix Mr to be arbitrary. 3 phases in Mr can be removed, leaving a total of 9 parameters 
which determine both the low energy neutrino masses and mixings as well as the baryon 
asymmetry. It then becomes apparent that T]b is calculable in terms of the neutrino observ- 
ables. There have been other attempts in the literature to relate leptogenesis with low energy 
observables jHJ E]. Such attempts often make additional assumptions such as Mr, = M up 
(which may not be fully realistic), or specific textures for lepton mass matrices. 

While a lot has been learned from experiments about the light neutrino masses and 
mixings, a lot remains to be learned. Our analysis shows that cosmology puts significant 
restrictions on the light neutrino parameters. Successful baryogenesis requires within our 
model that three conditions be satisfied: tan 2 9u — mi/m 2 , (3 ~ a + tt/2 and #13 = (0.01 — 
0.07). Here #12 and #13 are elements of the neutrino mixing matrix, m ; are the light neutrino 



mass eigenvalues and a, (5 are the Majorana phases entering in the amplitude for neutrinoless 
double beta decay. Future neutrino experiments will be able to either confirm or refute these 
predictions. 

The rest of the paper is organized as follows. In Sec. 2 we review briefly the minimal 
left-right symmetric model. In Sec. 3 we analyze leptogenesis in this model. Here we derive 
constraints imposed on the model from the requirement of successful leptogenesis. In Sec. 4 
we calculate the lepton asymmetry parameter e% generated in the model in Ni decay. Sec. 5 
summarizes the relevant Boltzmann equations needed for computing the baryon asymmetry 
parameter. Sec. 6 provides our numerical results for t)b- Finally, in Sec. 7 we conclude. 

2 Brief review of the minimal left-right symmetric model 

Let us briefly review the basic structure of the minimal SUSY left-right symmetric model 
developed in Ref. [7j. The gauge group of the model is SU(3) C x SU(2) L x SU{2) R x U(1) B -l- 
The quarks and leptons are assigned to the gauge group as follows. Left-handed quarks and 
leptons (Q,L) transform as doublets of SU(2)l [(5(3,2,1,1/3) and L(l, 2, 1, — 1)], while 
the right-handed ones (Q c , L c ) are doublets of SU(2) R [<5 C (3*, 1, 2, -1/3) and L c (l, 1, 2, 1)]. 
The Dirac masses of fermions arise through their Yukawa couplings to a Higgs bidoublet 
$(1,2,2,0). The SU(2) R x U(1) B -l symmetry is broken to U(1) Y by the VEV (v R ) of a 
B — L = —2 triplet scalar field A c (l, 1, 3, —2). This triplet is accompanied by a left-handed 
triplet A(l, 3, 1, 2) (along with A and A c fields, their conjugates to cancel anomalies). These 
fields also couple to the leptons and are responsible for inducing large Majorana masses for 
the v R . An alternative to these triplet Higgs fields is to use B — L = ±1 doublets ^(1, 2, 1, —1) 
and x c (l, 1,2,1), along with their conjugates x an d X c - m this case non-renormalizable 
operators will have to be invoked to generate large v R Majorana masses. For definiteness we 
shall adopt the triplet option, although our formalism allows for the addition of any number 
of doublet Higgs fields as well. The superpotential invariant under the gauge symmetry 
involving the quark and lepton fields is 

W = Y q Q T T 2 §T 2 Q c + YjL t t 2 $t 2 L c + (fL T ir 2 AL + f c L cT zr 2 A c L c ) . (3) 
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Under left-right parity symmetry, Q <-> Q C *,L <->• L c *,$ <-> $t, A <-» A c *, along with 

W SU (2) L ^ w *su(2) R i w b-l <-> Wb-£ and ^ ^ As a consequence, Y 9 = Yj, Y z = Yj, 
and f = f * in Eq. (jSJ). 4 It has been shown in Ref. [7^ that the hermiticity of the Yukawa 
matrices (along with the parity constraints on the soft SUSY breaking parameters) helps to 
solve the supersymmetric CP problem that haunts the MSSM. 

Below vr, the effective theory is the MSSM with its H u and Higgs multiplets. 5 These 
are contained in the bidoublet $ of the SUSY left-right model, but in general they can 
also reside partially in other multiplets having identical quantum numbers under the MSSM 
symmetry (such as the \i X doublet Higgs fields alluded to earlier). Allowing for such 
a possibility, the superpotential of Eq. (j3J) leads to the relations for the MSSM Yukawa 
coupling matrices 

Y u = 7 Y d , Y t = -yY u n . (4) 

These relations have been called up-down unification [7j. Here, the first relation of Eq. 
Pj) implies 2jt ~ 7 tan/? = c where 7 is a parameter characterizing how much of H u and 
H(i of MSSM are in the bidoublet The case of H u d entirely in $ will correspond to 
7 = 1 and tan/? = mt/rrib- At first sight the first of the relations in Eq. (J3J) might appear 
phenomenologically disastrous since it leads to vanishing quark mixings and unacceptable 
quark mass ratios. It was shown in the first paper of Ref. [7j that including the one-loop 
diagrams involving the gluino and the chargino and allowing for a flavor structure for the 
soft SUSY breaking A terms, there exists a large range of parameters (though not the entire 
range possible in the usual MSSM) where correct quark mixings as well as masses can be 
obtained consistent with flavor changing constraints. 

It is the second of Eq. (0J) that concerns us here. This relation would lead to Md = cMi, 
with c ~ mt/rrib. One can thus go to a basis where the charged lepton and the Dirac neutrino 
mass matrices are simultaneously diagonal. The heavy Majorana mass matrix Mr = ¥vr will 
then be a generic complex symmetric matrix. After removing three phases in Mr by field 
4 We do not explicitly use these relations. 

5 The right-handed gauge bosons have masses of order vr ~ 10 14 GeV and thus play no significant role 
in cosmology at T ~ Mi <C vr. 
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redefinitions, we are left with 9 parameters (6 magnitudes and 3 phases) which determine 
the light neutrino spectrum as well as the heavy neutrino spectrum. This in turns fixes the 
lepton asymmetry. The consequences of such a constrained system for leptogenesis will be 
analyzed in the next section. 

In principle the A(l,3,l,+2) Higgs field can also acquire a small VEV < 0(eV) [TT] . 
In this case the seesaw formula would be modified, as will the calculation of the lepton 
asymmetry IT] . We will assume such type II seesaw contributions proportional to (A) are 
zero in our analysis. This is consistent with the models of Ref. [7]. Leptogenesis in the 
context of more general left-right symmetric models has been analyzed in Ref. ^2] • 

3 Leptogenesis in left-right symmetric framework 

The SU(2) R x U(1) B -l symmetry is broken down to U(l) Y by the VEV (A c ) = v R ~ 
10 14 GeV. At least some of the right-handed neutrinos have masses below v r. We thus focus 
on the neutrino Yukawa coupling in the context of MSSM. The SU(2)l x U(1)y invariant 
Yukawa interactions are contained in the MSSM superpotential 



where / stands for the left-handed lepton doublet, and (e c , u c ) denote the conjugates of 
the right-handed charged lepton and the right-handed neutrino fields respectively. H u , H d 
are the MSSM Higgs fields with VEVs v u , v^. Mi = Y e va, M D = Y v d v u and M R are 
respectively the charged lepton, the Dirac neutrino, and the Majorana r.h.n mass matrices. 
Then one can generate light neutrino masses by the seesaw mechanism [2] 



There is mixing among generations in both M R and Md, the light neutrino mixing angles 
will depend on both of these mixings. Within the SM or MSSM where is an arbitrary 
matrix, the structure of the right-handed neutrino mass matrix can not be fully determined 
even if the light matrix M v were to be completely known from experiments. As noted in 



W = lH d Y t e c + IH u Y u d v c + - v' 




(5) 



M v = -M D M^M D T . 



(6) 



6 



Sec. 2, in the minimal version of the left-right symmetric model one has 



M. 



D 



cMi = c diag{m e ,m v ,,m T ) 



(7) 



where c ~ Here we have already gone to a basis where the charged lepton mass matrix 
is diagonalized. In the three family scenario, the relations between the flavor eigenstates 
(u e , Up, v T ) and the mass eigenstates (ui, z/ 2 , z/ 3 ) can be expressed in terms of observables as 



M v = U*Mf a9 U\ 



(8) 



where M^ iag = diag(mi,m2,m 3 ), with rrii being the light neutrinos masses and U being the 
3x3 mixing matrix which we write as U = Upmns-P- We parameterize Upmns [12] as 

(u el u e2 u e3 \ 



u 



PMNS 



\u Tl u t2 U t3 J 
( C12C13 

— S 12C23 ~~ c 12Sl3S23e l 



S12C13 
C12C23 ~~ s 12 s 13 s 23 el 



C13S23 



\ Sl2 s 23 — Ci2Si3C23e M — C12S23 — Si2S 13 C23e* <5 C13C23 J 

where c^- = cos%, Sij = sin dij and S is the Dirac CP violating phase which appears in 
neutrino oscillations. The matrix P contains two Majorana phases unobservable in neutrino 
oscillation, but relevant to neutrinoless double beta decay 

/ e ia \ 

P= e lf3 . (10) 

\o 1) 

Combining Eq. (jHJ) with the seesaw formula of Eq. © and solving for the right-handed 
neutrino mass matrix we find 



(9) 



M 



R 



c 2 M l M~ 1 M l 

/me 0\ 



c 2 ml 



mi 



0^0 

V \] 



(\ \ 



UpMNsP 2 



2* 







u 



PMNS 



\0 m 

\ U13 / 



/me 0\ 
0^0 

V I) 



.(ID 



This enables us to establish a link between high scale parameters and low scale observables. 



We define a small expansion parameter e as 



e = 



hi,, 

— ^ ~ 0.059, 



in terms of which we have 



m e = a e e 3 m. 



rrii 

— = a 13 e, 



! is = t 13 e, e. 



23 = 




(12) 



Here a e , a i3 , t i3 and t 23 are < O(l) parameters with a e = 1.400. These expansions follow 
from low energy data assuming the picture of hierarchical neutrino masses. 

We find that the requirement of generating adequate baryon asymmetry places significant 
constraints on the neutrino mixing parameters. Specifically, the following expansions 



where ai 2 and b are < 0(1) parameters are required. To see this, we note that the CP 
asymmetry parameter e\ generated in the decay of N\ is too small, of order e\ ~ ~ 2 x 10~ 9 
if a\2 or b are much greater than 1. This is because the heavy neutrino masses would be 
strongly hierarchical in this case, M\ : M2 : M3 ~ e 6 : e 2 : 1. This can be altered to a weak 
hierarchy M\ : M 2 : M 3 ~ e 4 : e 2 : 1 by observing that the elements of the 2-3 block of Mr 
of Eq. ()11|) are all proportional to {^e 2l/3 cos 2 # 12 + e 2ia sin 2 # 12 } and by demanding this 
quantity to be of order e. Eq. (jT3j) is just this condition. S\ ~ |^ ~ 10 -6 in this case, which 
can lead to acceptable baryon asymmetry, as we show. 

An immediate consequence of Eq. (fTB^) is that neutrinoless double beta decay is sup- 
pressed in the model. The effective mass relevant for this decay is found to be 



This is of the order m 3 e 2 ~ 10~ 4 eV, which would be difficult to measure. This amplitude 
is small because of a cancelation between the leading contributions proportional to ni\ and 
m 2 (see Eq. (fT3"jO. 

In terms of these expansions, the r.h.n mass matrix becomes 




(13) 




(14) 



Mr 



/A n e 5 A 12 e 3 A 13 e 2 \ 

A 12 e 3 A 22 e 2 A 23 e 
\A 13 e 2 A 23 e A 33 J 



(15) 
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where 



A 

A 

A 

A 
A 



11 



12 — 



13 — 



-22 — 



M ea 2 e e 2ia cos2# 12 

COS 2 6*12 

M ea e e 2ia tan6 12 



23 — 



V2 

M ea e e 2ia tand 12 

^ {a 13 - a 12 e 2ia cos 2 i2 - 2ibe 2ia sin 2 12 + 2e* (2a+<5) t 13 tan# 12 } 
(ai 3 + a l2 e 2ta cos 2 12 + 2^e 2ia sin 2 9 12 } 



A 33 = -M Q ee 



2ia 



2 2 

Here we defined M Q = This 

mi 



J, l5+ fl , , . 2/3 ai 3 e 2m ai 2 cos 2 gi 2 ) 
<^ t 13 e w tan0 12 + 2&sm z 12 1 > . 



2 

'-. This hierarchical mass matrix is diagonalized by a series 

nil ° J 

rotations Ui, U 2 and U 3 such that; 

[{Mil \ 

{KU^U^Mr^KU^U: 



(16) 
ies of 



'i) T = 



|M 2 | 
V \M 3 \J 



(17) 



where .K" = diag(ki, k 2 , k 3 ) with ki = e~ l<t> ^ 2 being phase factors which make each r.h.n 
masses M; real, M; = |M,|e^. V = (KU 3 U 2 U 1 ) T is the matrix that diagonalizes M fl . The 
unitary matrix U\ is given by 

/ 1 Tt' 2 \ 



Ui = 



1 

1 ) 



(18) 



v ^33 ' 

Similarly, U 2 and U 3 are unitary matrices with off-diagonal entries given by 



(^2)23 = -4^e , (^3)12 = 
^33 



4 ^13^23 ^ , 

A12 — ) e 



The mass eigenvalues are found to be 



4 _ ^23 

^ 22 AiF 



(19) 



Mi = M Q k 2 e 5 (2a 13 a 2 e e 2ta sin 2 6 12 ) 

x (2t 2 3 e 2l(a+<5) sin 2 9 12 + (a i2 + 2«6 + (a 12 - 2ib) cos 29 12 )a 13 cos 2 0i 2 ) _1 
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M 2 = M kje 3 e 2ia (a 13 (a 12 + lib + (a 12 - lib) cos 2B 12 ) + 2tj 3 e t{S+a) tan 2 12 ) 
x (-ai 3 + «&e 2ia + e 2ia (ai 2 cos 2 ^2 -^cos20 12 ) + 2e t<5 t 1 3tan^ 12 ) _1 



M 3 = — p- (a 1 3-^e 2 ^-e 2M (ai 2 cos 2 ^i 2 -z6cos2e 12 + 2t 13 e i5 tan^ 2 )). (20) 
We use these results in the next section to determine e\. 



4 CP violation and lepton asymmetry 

Now that we have developed our framework, we can turn attention to the evaluation of 
the CP asymmetry E\ generated in the decay of the lightest r.h.n N\. This arises from the 
interference between the tree-level and one-loop level decay amplitudes. 6 In a basis where 
the r.h.n mass matrix is diagonal and real, the asymmetry in the decay of iVj is given by [T5] 



' ^lm[(MlM D ) 12 



(21) 



3=2,3 

where f(x) and g(x) represent the contributions from vertex and self energy corrections 
respectively. For the case of the non-supersymmetric standard model with right-handed 
neutrinos, these functions are given by ^5] 

-1 + (X + I)\n[ 1 + -]] , Snon-SUSY(X) = , (22) 



/non-SUSY(^) — \/X 



X 

while for the case of MSSM plus right-handed neutrinos, they are given by 

(1 \ 2 y/x 
1 + - ) , S'susy^) = . (23) 

X J x — 1 

Here v is the SM Higgs doublet VEV, v ~ 174 GeV. For the case of MSSM, v in Eq. (j2U) is 

replaced by v sin (3. Hereafter, for definiteness in the numerical evaluation of the Boltzmann 

equations, we assume the SM scenario. However, our result should be approximately valid 

for the MSSM case as well. 7 Assuming a mass hierarchy M\ -C M 2 < M 3 in the right-handed 

6 We will assume Mi -C M2 < M3. In this case, even if the heavier right-handed neutrinos N2 and N3 

produce lepton asymmetry, it is usually erased before the decay of iVi. 

7 The function / + g in MSSM is twice as big compared to the SM. However this is compensated by the 

factor j- that appears in t]b which in MSSM is half of the SM value. 
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neutrino sector i.e., (x ^> 1), which is realized in our model, see Eq. ([To]) , the above formula 
is simplified to the following one: 

£l = " ifi 2,Mi M , E M{M\,M D )l k \ ^. (24) 

Ei depends on the (1,1), (1,2) and (1,3) entries of M^ D M D . These quantities can be related 
to the light neutrino mass and mixing parameters measurable in low energy experiments. In 
the basis where Mr is diagonal, these elements are 

(i\4M D ) n = (cm T ) 2 (V 31 V* 1 + V 21 V* 1 e 2 + a 2 e V 11 V^) 
(MlM D ) 12 = {cm T ) 2 {V zl V; 2 + V 2l V* 2 e 2 + a 2 e V ll V* 2 e 6 ) 
(i\4M D ) 13 = {cm T ) 2 (V 31 V* 3 + V 21 V* 3 e 2 + a 2 e V 11 V 1 * 3 e 6 ), (25) 

where V = KU 3 U 2 Ui is the unitary matrix diagonalizing Mr. Straightforward calculations 
give, to leading order in e, 

(M],M D ) n = 8a 2 c 2 m 2 e 4 cos 2 6 12 sin 2 9 12 (a 2 13 + t 2 13 tan 2 9 12 ) 

x 1 / {8^3 sin 4 6 12 + 32a 13 t 2 13 b cos 2 6 12 sin 4 12 sin 2(a + 5) 

+ ai3 cos 4 #i 2 [4ai 3 (a 2 2 — b 2 ) cos 26 12 + ai 3 (a 2 2 + 46 2 )(3 + cos4#i 2 ) 

+ 16ai 2 t 2 3 sin 2 0i2Cos2(a + 5)]} (26) 



12 

2 



(M j D M D )\ 2 = 2a 2 c 4 m 4 e 6 tan 2 o 12e ~<4>i-4>2) e -M2a+5) _ t y cog 2 ^ i2 _ ^ gin 2 # 

(2ai 3 e l(2Q+5) - (a i2 + 2«6) e - i(5 ) + 4(a 2 3 + t 2 3 ) + £ 13 sin40 12 e-* 5 (ai 2 - 2ib)} 
x 1/ { [ibe lS - a 13 e- l{2a+S) + a 12 e~ lS cos 2 6 12 - %be' l& cos 26 12 + 2* i3 tan #i 2 ] 2 
x [3ai 2 ai 3 - 2ia l3 b + At 2 13 e- 2l( - a+s) + 4 cos 26 12 (a 12 a 13 - t\ 3 e- 2l{a+S) ) + 

ai 3 (ai 2 + 2^)cos4#i 2 ] 2 } (27) 



2a 2 c 4 m 4 e 4 sin 2 9 l2 e-^-^\a l3 cos6 l2 + e' l{2a+5 h 13 sin6 12 ) 2 
1/ {ai 3 cos 2 M°i2 - 2?& + (a i2 + 2^6) cos 20 12 ) + 2t\ 3 sin 2 #i 2 e- 2 * (Q+5) } 2 (28) 

These analytical expressions have been checked numerically. In Figure (1) we have plot- 
ted |ei| as function of 6*i 3 for fixed values of other observables. The solid line in Fig (1) 
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(MlM D ) 2 13 = 

x 



0.05 0.1 0.15 0.2 

Figure 1: Plots for CP asymmetry parameter e% using analytical (dotted) and numerical 
(solid) results as a function of the neutrino oscillation angle #13. The input parameters used 
are a 12 = 1, b = 1, Am| = 2.5 x l(T 5 eV 2 , Am 2 a = 5.54 x 10~ 3 eV 2 and {5, a} = {n / 4, tt / '4} . 
Our model requires \si\ > 1.3 x 10~ 7 to successfully generate an adequate number for the 
BA. This criterium happens to be satisfied only in the region for which 0.01 < #13 < 0.07, 
this interval is not too sensitive to variations in the input parameters. 

which corresponds to the exact numerical evaluation agrees very well with the dashed line 
corresponding to the analytical expressions. 

From Figure (1), it is apparent that #13 is constrained in the model from cosmology. If 
£1 < 1.3 x 10~ 7 , the induced baryon asymmetry would be too small to explain observations. 
As can be seen from Figure (1), #13 should lie in the range 0.01 — 0.07 for an acceptable value 
of £1. This result does not change very much with variations in the other input parameters. 

Electroweak sphaleron processes [BJ will convert the induced lepton asymmetry to baryon 
asymmetry. The ratio of baryon asymmetry to entropy Yg is related to the lepton asymmetry 
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through the relation [To] : 



CY, 



C 



Y L 



(29) 



B—L — 



C-l 



where C = 22N +13N ' = 3 and N v = 1, 2 in the case of the SM and MSSM respectively. 

In either case C ~ §. In Eq. (123), Vb = ^ with s = 7.04 n T 

There has been considerable interest in obtaining approximate analytical expression for 

baryon asymmetry |17| ITS] . In order to estimate this, the dilution factor, often referred 

to as the efficiency factor k that takes into account the washout processes (inverse decays 

and lepton number violating scattering) has to be known. As an example, k = (2 ± 1) X 
/ y\i.i±o.i 

/ o.okv j j-^g b een suggested in Ref. [T7j from which r/ B ~ 0.96 x lO -2 ^^ has 



been calculated. In our work we solve the coupled Boltzmann equations numerically to 
estimate the baryon asymmetry without referring to the efficiency factor. 

5 Boltzmann equations 

In this section we set up the Boltzmann equations for computing the baryon asymmetry 
rjB generated through the out of equilibrium decay of N\. In our model the right-handed 
neutrino masses are not independent of the CP asymmetry parameter e±. Therefore a self 
consistent analysis within the model is required. 

In the early universe, at temperature of order N\ mass, the main thermal processes which 
enter in the production of the lepton asymmetry are the decay of the lightest r.h. neutrino, 8 
its inverse decay, and the lepton number violation scattering, AL = 1 Higgs exchange plus 
AL = 2 r.h.n exchange ^H]- The production of the lepton asymmetry via the decay of 
the r.h.n is an out-of-equilibrium process which is most efficiently treated by means of the 
Boltzmann equations (BE). 

The first BE which describes the evolution of the abundance of the r.h. neutrino and 
8 In our analysis we stick to the case where the asymmetry is due only to the decay of the lightest r.h. 
neutrino N±. 
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which corresponds to the source of the asymmetry is given by 9 

%=^GM^ + ^ (30) 

where z = Here s(z) is the entropy density and j D , j s are the interaction rates for the 
decay and AL = 1 scattering contributions, respectively. 

The second BE relevant to the lepton asymmetry is given by 

dz ' s(z)H{M 1 ) [ £i1d ^ \Y^ q J + lw Y e L q \ ' 1 ' 

where e\ is the CP violation parameter given by Eq. (|2T|) and j w is the washout factor 
which is responsible for damping of the produced asymmetry, see Eq. (|49|) below. In Eqs. 
()30|) and 1)31)1. F/ 9 is the equilibrium number density of a particle species z, which has a mass 
mj, given by 

Y*(z) = ^*(^X*kJ^\ , (32) 



4vr 4 ^ V^i7 ~ \M X 
where g>j is the particle internal degree of freedom (g N = 2, g e = 4). At temperatures far 
above the electroweak scale one has g* ~ 106.75 in the standard model, and g* ~ 228.75 
in MSSM. H, the Hubble parameter evaluated at z — 1, and the entropy density, are 
given by 

^V^rs; • •<*> = - sr 1 -?-- < 33) 

where M P = 1.22 x 10 19 GeV. We also have 

T Sj =2 7 W+47f. (34) 
The decay reaction density 7 D has the following expression: 

where K n (z) are the modified Bessel functions. ITV of the r.h.n Nj is the tree level total 
decay rate defined as 

T Nj = ^M, (36) 



9 In this section we follow the notation of the first paper of Ref. )S] to which we refer the reader for further 
details. 
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where 



ei(rr\ diTrrii (m 



2vr 2 



T 



(37) 



We used the definition A = M D /v. We define the reaction density 7M of any process 
a + b c + d by 

,(,) _ M? 1 



7 V 



647T 4 2 j (Mg+M b ) 



dx a^\x)y/x K\ i\fxz). 



(3* 



where <7^(x) are the reduced cross sections for the different processes which contribute to 
the Boltzmann equations. For the AL = 1 processes involving the quarks, we have 



a 



(i) 



a 



(2) 



3a u E (Ktj^aj) 
a=l 
3 

3««E ( A «j Ac *i) 



where 



a=l 



a. 



00 CL a 



x 



X 



x — 2a,- + 2dh dj — 2ah f x — dj + a,h 
H — In 1 ' 



X — dj ■ + dh X — dj 



a h 



Alt 



Mi 



M l 



(39) 
,(40) 

(41) 



fi is the infrared cutoff which we set to 800 GeV fHl [2D]- For the AL = 2 r.h.n exchange 
processes, we have 

33 33 



a 



(i) 

N 



a 



(2) 
N 



E E ( A ^) ( A ^) 4- + E E ^ e ( A «n^) (a^) 2$ (42) 

a=l j=l a=l n<j,j=X 

3 3 3 3 

E E ( A ai A °i) ( A aj A "i) + E E Re ( A «n A "i) ^ ( 43 ) 

o=l j=l a=l n<j,j=l 



where 



A?) = - 



03 



B 



n.) 



1 

2^ 
1 

1~TX 



1 + — + — - 1 + 

Dj 2D] x 1 



x + a, 



In 



X + dj 



x 



+ 



a , 



x + dj x + 2dj 



In 



^Z/ I CL j 



(1) _ V^n C 



27T 



X 



+ + DjD n 
1 



+ 1 



3lj ' \ CXiy^ CL j \-) y 



111 



(1 + -JI ^Hn 

x / \dj — d n Dj 



X I 



a, 



5 



(2) 

7iJ 



2tx \ x + d n + dj 



In 



(x + dj)(x + a n ) 



CL j Q"yi 



2 , / d n (x + dj 
H In' v 



Oj^l CXi r-j \^ j I CXiiy^ J 



(44) 
(45) 
(46) 



(47) 
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and 

_ (x - dj) 2 + djCj _ ^ {KjKjKjKj) (A , 



x — a,- - - 647r 2 

J a=l 



Finally, ^ w that accounts for the washout processes in the Boltzmann equations is 



t ( k + SW," + 2 ^> - ^ ) + < + - (49) 



1w 

3=1 



9 id j Y 
Z 1 N. 



Here, we emphasize the so-called RIS (real intermediate states) in the AL = 2 interactions 
which have to be carefully subtracted to avoid double counting in the Boltzmann equations. 
This corresponds to the term — |7 D in Eq. P9*|). For more details see Refs. [T7||2I] and the 
first paper of Ref. 



6 Results and discussion 

We are now ready to present our numerical results. First we make several important remarks. 
Even though our model is supersymmetric, we have considered in our BE analysis only the 
SM particle interactions. This is a good approximation (see footnote 7). The authors in Ref. 
[2*U] have demonstrated that SUSY interactions do not significantly change the final baryon 
asymmetry. Furthermore, we have not included in our analysis the effects of renormalization 
group on the running masses and couplings. The first paper of Ref. [TH] has studied these 
effects. This paper has also included finite temperature effects and AL = 1 scattering 
processes involving SM gauge bosons, which we have ignored in our analysis. This should 
be a good approximation since it is believed that these effects are significant in the weak 
washout regime and our model parameters seem to favor the strong washout regime with 
rhi = {MdMbIh ~ o.l eV. Scattering processes involving gange bosons have also been 
studied in Ref. f2T] in the context of resonant leptogenesis where they have been shown to 
be significant. 

Our next step is to put this model to the test and check its predictions. In order to 
compute the value of the baryon asymmetry we proceed to numerically solve the Boltzmann 
equations. We scan the parameter space corresponding to the parameters a 12 , b, the oscilla- 
tion angle #13, the CP phase 5 and the Majorana phase a. In order to automatically satisfy 
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the oscillation data, we input the following light neutrino parameters: 

Am 2 Q = 2 - 5 x 10" 5 eV 2 , Am 2 = 5.54 x 10" 3 eV 2 , sin# 12 = 0.52. (50) 

Using hierarchical spectrum, we see that the masses mi, m 2 and m 3 are fixed. On the 
other hand we consider maximal mixing in the 2-3 sector of the leptonic mixing matrix, i.e 
$23 — f +t 2 3€ with t 23 being zero ( t 23 ~ C(l) has minimal impact on t]b)- The CP phase 5 
and the Majorana phase a are allowed to vary in the intervals [0, 2ir] and [0, tt] respectively. 
We remind the reader that the second Majorana phase (3 is related to a through (3 ~ a+^+be. 
6 13 will be allowed to vary in the interval [0; 0.2] as it is bounded from above by reactor 
neutrino experiments. 

In Figure (2), for a given set of input parameters, we illustrate the different thermally 
averaged reaction rates T„ = -^£- contributing to BE as a function of z = ^f. 

All rates at z — 1 fulfill the out of equilibrium condition (i.e. Tx J$ H(z = 1)), and so the 
expected washout effect due to the AL = 2 processes will be small. The parameters chosen 
for this illustration are: 5 = ir/2, a = tt/2, ai 2 = 0.01, b = 0.9, cm T = m t (j^j — 135 GeV 
and 6 n = 0.02. Eq. (jHOD fixes the light neutrino masses to be: m x = 0.00271292 eV, 
m 2 = 0.00688186 eV and m 3 = 0.0380442 eV. For this choice we obtain | e x |~ 2 x 10" 7 . 
The calculated r.h.n masses in this case are 

Mi = 9 x 10 9 GeV, M 2 = 8.7 x 10 11 GeV, M 3 = 2.6 x 10 14 GeV. (51) 

The mass of the lightest r.h.n is consistent with lower bound derived in Ref. JH] 5 M\ > 
2.4 x 10 9 GeV, for hierarchical neutrino masses assuming that one starts with zero N\ initial 
abundance (which is what we assumed in our calculation). This mass is also in accordance 
with the upper bound found in Ref. [22] following a model independent study of the CP 
asymmetry, and the bound derived in Ref. ^7j based on the estimation of z/r production 
and the study of the asymmetry washout. 

Figure (3) represents the solution of the BE, N\ abundance and the baryon asymmetry 
both as functions of z for the same set of parameters mentioned above. The final baryon 
asymmetry, in terms of the baryon to photon ratio, is (see dark, solid curve in Fig. (3) for 
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z» 1) 

r] B ~ 6.03 x 1(T 10 . (52) 

This number is inside the observational range of Eq. (0). Our codes were tested to reproduce 
the results in the first paper of Ref. |Sj before being applied to this model. 

At this point we should mention some potential difficulties with gravitino decays in the 
model. If supersymmetry breaking is mediated by conventional supergravity, it is natural 
to expect the mass of the gravitino to be of order 1 TeV. Such a gravitino, with its Planck 
scale suppressed interactions, would decay into MSSM particles with a lifetime of order 1 
second. The decay products would upset the success of standard big bang nucleosynthesis. 
In an inflationary scenario, demanding consistency with light element abundance puts a 
lower limit on the reheating temperature Tr. For gravitino mass in the range 300 GeV to 3 
TeV, a limit Tr > 6 x 10 6 GeV has been derived |23]. Now, the decays of N\ that generates 
lepton asymmetry should occur after reheating, since any asymmetry produced prior to that 
will be diluted by inflation. Thus Mi < Tr is required, which for gravitino mass in the range 
300 GeV to 3 TeV is in conflict with the predictions of Eq. (|5ip. 

There are several ways around this problem, (i) In gauge mediated SUSY breaking 
scenario the gravitino is the lightest SUSY particle with mass in the range 10 -4 eV— 100 GeV. 
For nig < 100 MeV, there are no cosmological or astrophysical problems. In such a scenario 
the axion can serve as the dark matter, (ii) In anomaly mediated SUSY breaking scenario, the 
gravitino mass is enhanced by a loop factor compared to the squark masses and is naturally 
of order 100 TeV. Such a gravitino would decay with a shorter lifetime without affecting big 
bang nucleosynthesis. The gaugino is a natural dark matter candidate in this case, (iii) The 
gravitino itself can be the LSP and dark matter with a mass of order 100 GeV, in which case 
it does not decay [21]. Other solutions include changing the dynamics of the leptogenesis 
process by invoking (iii) non-thermal leptogenesis j2H] > (iv) resonant leptogenesis |2T) 125]. 
or (v) soft leptogenesis [27] . The results presented here are compatible with any one of the 
first three alternatives. 
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0.01 0.1 1 10 100 

z = Mi/T 



Figure 2: Various thermally averaged reaction rates Tx contributing to BE normalized to the 
expansion rate of the Universe H (z = 1). The straight greyed line represents H(z)/H(z = 1), 
the dashed is for Tu 1 /H(z = 1), the dotted-dashed line represents Tal=i/ H(z = 1) processes 
and the red curve represents Tal=2/H(z = 1). 



19 




0.01 0.1 1 10 100 

z = m Nl /T 



Figure 3: Evolution of Y Nl (solid blue), (dot-dash) and the baryon asymmetry r] B (dark 
solid line) in terms of z in the model. The estimated value for the baryon asymmetry is 
i]b — 6.03 x 10 -10 , with Y™ 1 = and assuming no pre-existing B — L asymmetry. 
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7 Conclusion 

An attractive feature of the seesaw mechanism is that it can explain the origin of small 
neutrino masses and at the same time account for the observed baryon asymmetry in the 
universe by the out of equilibrium decay of the super-heavy right handed neutrinos. It is 
then very tempting to seek a link between the baryon asymmetry parameter t]b induced 
at high temperature and neutrino mass and mixing parameters observable in low energy 
experiments. No quantitative connection can be found between them in the SM. There have 
been several attempts in the literature [HI EJ I2E] to establish a relationship between the two. 
In this paper we have addressed this question in the context of a class of minimal left-right 
symmetric models. 

In the models under consideration the minimality of the Higgs sector implies that Mi and 
Md (charged lepton and Dirac neutrino mass matrices) are proportional. As a result, the 
entire seesaw sector (including the heavy right-handed neutrinos and the light neutrinos) has 
only 9 parameters. This is the same number as low energy neutrino observables (3 masses, 
3 mixing angles and 3 phases). As a result we are able to link the baryon asymmetry 
of the universe to low energy neutrino observables. This feature is unlike the SM seesaw 
which has too many arbitrary parameters. Our numerical solution to the coupled Boltzmann 
equations shows that this constrained system with Mi cx Md leads to an acceptable baryon 
asymmetry. The requirement of an acceptable baryon asymmetry restricts some of the light 
neutrino observables. We find that tan 2 8 12 ~ mi/m 2 , 0.01 < #13 < 0.07 and (3 ~ a + 7r/2 
are needed for successful baryogenesis. Future neutrino oscillation experiments can directly 
probe into the dynamics of the universe in its early stages. 
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